[LONGITUDINAL]	
input={K, KDE, KPQ, KQPP, LAMBDAP, GAMMA, DELTAQP}

PK:
depot(target=C)

EQUATION:
C_0 = 0
PT_0 = 5
Q_0 = 40
QP_0 = 0
Q_0 = 40
PSTAR = PT+Q+QP
ddt_C = -KDE*C
ddt_PT = LAMBDAP*PT*(1-PSTAR/K) + KQPP*QP - KPQ*PT - GAMMA*KDE*PT*C
ddt_Q = KPQ*PT - GAMMA*KDE*Q*C
ddt_QP = GAMMA*KDE*Q*C - KQPP*QP - DELTAQP*QP

DEFINITION:
PSTAR_obs = { distribution=normal, prediction=PSTAR,  errorModel=constant(a)}


[INDIVIDUAL]
input={pop_LAMBDAP,pop_KQPP,omega_KPQ,pop_KDE,omega_KDE,pop_DELTAQP,omega_DELTAQP,omega_KQPP,
       pop_GAMMA,pop_KPQ,omega_GAMMA,omega_LAMBDAP}

DEFINITION:
LAMBDAP = { distribution = lognormal, reference = pop_LAMBDAP, sd=omega_LAMBDAP}
KQPP =    { distribution = lognormal, reference = pop_KQPP, sd=omega_KQPP}
KPQ =     { distribution = lognormal, reference = pop_KPQ, sd=omega_KPQ}
KDE =     { distribution = lognormal, reference = pop_KDE, sd=omega_KDE}
GAMMA =   { distribution = lognormal, reference = pop_GAMMA, sd=omega_GAMMA}
DELTAQP = { distribution = lognormal, reference = pop_DELTAQP, sd=omega_DELTAQP}




